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Most cosmological structures in the universe spin. Although 
structures in the universe form on a wide variety of scales 
from small dwarf galaxies to large super clusters, the gen- 
eration of angular momentum across these scales is poorly 
understood. We have investigated the possibility that fila- 
ments of galaxies - cylindrical tendrils of matter hundreds 
of millions of light-years across, are themselves spinning. By 
stacking thousands of filaments together and examining the 
velocity of galaxies perpendicular to the filament’s axis (via 
their red and blue shift), we have found that these objects 
too display motion consistent with rotation making them 
the largest objects known to have angular momentum. The 
strength of the rotation signal is directly dependent on the 
viewing angle and the dynamical state of the filament. Just 
as it is easiest to measure rotation in a spinning disk galaxy 
viewed edge on, so too is filament rotation clearly detected 
under similar geometric alignment. Furthermore, the mass 
of the haloes that sit at either end of the filaments also in- 
creases the spin speed. The more massive the haloes, the 
more rotation is detected. These results signify that angular 
momentum can be generated on unprecedented scales. 


Introduction 


How angular momentum is generated in a cosmological con- 
text is one of the key unsolved problems of cosmology. In 
the standard model of structure formation, small overdensities 
present in the early universe grow via gravitational instability as 
matter flows from under to overdense regions. Such a potential 
flow is irrotational or curl-free: there is no primordial rotation 
in the early universe and angular momentum must be generated 
as structures form. 


Tidal torque theory Hl provides one explanation - the mis- 
alignment of the inertia tensor of a gravitationally collapsing re- 
gion of space with the tidal (shear) field can give rise to torques 
which spin up the collapsing material [BE] Such an explanation 


is valid only in the linear regime namely in limit where density 
perturbations are small with respect to the mean and where flows 
are laminar. As a collapsing region reaches turn-around tidal 
torques cease to be effective and the final angular momentum 
of a collapsed region is far from what tidal torque theory would 
predict ELN Although one recent study has demonstrated that 
galaxy spin direction (i.e. clockwise versus counter clockwise) 
can be predicted from initial conditions, revealing a critical clue 
to the non-linear acquisition of angular momentum, our under- 
standing of spin magnitude, direction and history remains in its 
infancy. Regions that are still in the linear or quasi-linear phase 
of collapse could provide a better stage for the application of 
tidal torque theory. 


Cosmic filaments [2 being quasi-linear extended topograph- 
ical features of the galaxy distribution, provide such an environ- 
ment. Yet, owing to the challenges in characterizing and iden- 
tifying such objects, potential rotation on the scales of cosmic 
filaments has been discussed! but never measured until now. 


It is known that the cosmic web in general and filaments, in 
particular, are intimately connected with galaxy formation and 
evolution LE] They also have a strong effect on ee 
16l often regulating the direction of how galaxies and 
their dark matter halos rotate 52/41) However, it is not known 
whether the current understanding of structure formation pre- 
dicts that filaments themselves, being uncollapsed quasi-linear 
objects should spin. A recent study A eed on the ArXiv 
while this draft was being finalized) “=! examined the velocity 
field around galactic filaments defined by halo pairs in a large 
N-body simulation and found a statistically significant rotation 
signal. This is an intriguing finding and, although filaments and 
their rotation speed are defined differently, the current work is 
in which the observed galaxy distribution is examined in a bid 
to find possible filament rotation was partly motivated by the 
theoretical suggestion that filaments may spin a 


Results 


After segmenting the galaxy distribution into filaments us- 
ing a marked point process known as the Bisous model) each 
filament can be approximated by rectangle on the sky and thus 
the galaxies within it may be divided into two regions on either 
side of the filament spine. The mean redshift difference Azap of 
galaxies between two regions are considered as a proxy for the 
line of sight velocity difference and hence for the filament spin 
signal. Any measured value of Azag needs to be assigned a sig- 
nificance based on a randomization procedure (explained in the 


methods section). Figure[I]shows the statistical significance of 
the measureed Azap as a function of Zms/Azap - a proxy for 
the dynamical “temperature” of the filament (see method sec- 
tion for more details). The number of galaxies in each region is 
denoted by color. Two salient points can be gleaned here. First, 
The more galaxies in a given filament, the more inconsistent the 
redshift difference Azap is with random. Second, (as expected) 
the colder the filament, the more inconsistent the redshift dif- 
ference is with random. This second point is a generalization - 
cold filaments with zms/Azap < 1 and few galaxies, can have 
redshift differences only weakly inconsistent with random ex- 
pectations. However one may note that as a trend, the colder 
the filament, the more significant the redshift difference is. In 
other words if Azap is considered a proxy for filament spin, we 
observe a spectrum of filaments from dynamically hot that are 
consistent with random to dynamically cold filaments that are 
completely inconsistent with random at the many sigma level. 
Note that even for dynamically hot filaments, there are a few 
that are highly inconsistent with random. The reader will note 
that the wide distribution of significance seen at a given value of 
Zrms /Azap, is also a reflection of the distribution of inclination 
angle made by the filament axis with the line of sight (see figure 
Al). 

The cumulative distribution of Azag for both the entire ob- 
served filament sample as well as the randomized trials, is 
shown in the three panels of Figure |2| The cumulative distri- 
bution is shown for all filaments (left), filaments whose axis is 
inclined by cos @ < 0.5 to the line of sight (middle) and dynam- 
ically cold filaments (Z;ms/Azap < 1) inclined by cos@ < 0.2 
(right). The reader will note that even when examining all fila- 
ments where inclination angle is completely ignored (and hence 
includes filaments viewed along their axis which will likely 
weaken the signal), the full distribution of Azag is inconsistent 
with randomization tests. This inconsistency increases when 
considering filaments with cos@ < 0.5 and cos@ < 0.2 that are 
cold. To quantify the statistical significance of the cumulative 
distribution one may simply measure, in units of the 10,000 ran- 
dom trial’s standard deviation, how far the measured signal is 
from the mean randomized signal. This is plotted in the upper 
panels of Figure [2] and shows that randomized trails are statis- 
tically inconsistent with the measured signal at very high con- 
fidence. In other words, shuffling the redshifts of galaxies in a 
filament is unlikely to produce redshift differences as great as 
that observed. 

The redshift difference (a proxy for the rotation signal) 
stacked across various (sub-) samples, is presented in Figure [3] 
The following conventions have been adopted. Region A (de- 
fined as the region with greater mean redshift) is plotted in the 
upper part of each plot, while region B is plotted in the lower 
part. The position of each galaxy along and perpendicular to 
the filament axis is shown on the x and y axis, respectively and, 
along the x axis, is normalized to the filament’s length. Each 
galaxy is colored by its redshift difference Az, with respect to 
the mean redshift of all galaxies in the filament according to 
the color bar at right. In the ideal situation where all galaxies 
exhibit circular or helical motion about the filament axis, such a 
plot would only have red points in the upper part and blue points 
in the lower part. The statistical significance of each (sub-) sam- 


ple is denoted by the significance, in units of o and indicated on 
top of each panel. 

In Figure [3] left column (i.e. panels, A, D, G) we show the 
stacked rotation signal for all filaments, filaments whose axis 
subtends an angle cos@ < 0.5 with the line of sight and fila- 
ments whose axis subtends an angle cos@ < 0.2 with the line 
of sight and which have zms/Azag < 1. Comparing panel (A) 
with (D) the reader will note what has been mentioned before 
namely that merely changing the inclination angle increases the 
signal. Panel (G) shows a very strong rotation signal - at 3.30 
- when considering dynamically cold filaments that are mostly 
perpendicular to the line of sight the rotation signal becomes 
very convincing. 

Since filaments are long tendrils of galaxies often connect- 
ing nodes of the cosmic web, the mass of the two halos closest to 
the filament’s two end points are examined for a possible corre- 
lation with the signal strength we measure. We thus perform the 
following procedure. We search the group catalogue #4! for the 
two groups that are closest to the filament’s two end points. The 
reader is referred to that paper for details. In sum, the galaxy 
groups are detected with a modified Friend-of-Friend (FoF) al- 
gorithm, where the FoF groups are divided into subgroups so 
that each subgroup can be considered as a virialized structure. 
The mass of each group is estimated using the virial theorem, 
the size of the group in the sky and the group velocity dispersion 
along the line of sight. By summing these two group masses we 
can ascribe a single mass abutting each filament. We call this 
mass the “filament end point mass”. 

The three samples of filaments can then be subdivided ac- 
cording to this characteristic. Figure[3|middle and right columns 
show the stacked signal for the 10 percent of the filaments 
with the smallest filament end point mass (ie panels B,E,H) and 
largest 10 percent filament end point mass (panels C, F, I) asso- 
ciated halo end points, respectively. The limits of the 10th and 
90th percentile group masses of panel H and I are ~ 10!°4 and 
~ 10143, in units of Mo, respectively. 

A clear dependence on filament end point mass is seen for 
all three samples in Figure[3| The smaller the filament end point 
mass, the weaker the rotation signal. The filaments with the 
largest end point mass, most perpendicular to the line of sight, 
and dynamically cold, show an outstanding 4.20 inconsistency 
with the null hypothesis of random origin. The fact that the mass 
of the clusters at the end of each filament directly affects the 
measured signal, is evidence that the signal we are measuring is 
a physical effect related to gravitational dynamics. 

A “rotation curve” can be constructed for the stacked fila- 
ment sample by converting the redshift difference into a velocity 
difference computed as c x Az;, where Az; is the redshift differ- 
ence between all the galaxies at given distance and the mean 
redshift of all galaxies in the filament. This is shown in Fig. 4] 
We adopt the typical convention that negative (positive) veloci- 
ties corresponds to approaching (receding) motion. Fig. [4]shows 
that rotation speed increases with distance, reaching a peak of 
around ~ 100km /s at a distance ~ 1 Mpc. It then decreases with 
distance, tending to zero at a distance of greater than ~ 2 Mpc. 
Incidentally this behavior motivates the choice of filament thick- 
ness in this study, and should be kept in mind in future studies 
of the dynamics of cosmic filaments. 


Summary & Discussion 


The results presented in this paper are consistent with the de- 
tection of a signal one would expect if filaments rotated. What is 
measured and presented here is the redshift difference between 
two regions on either side of a hypothesized spin axis that is 
coincident with the filament spine. The full distribution of this 
quantity (Figure[2) is inconsistent with random regardless of the 
viewing angle formed with the line of sight but does strengthen 
as more perpendicular filaments are examined. It also strength- 
ens when considering dynamically cold filaments - filaments 
whose galaxies have small rms values of their redshift. We em- 
phasize that this is not a trivial finding. There is no reason — be- 
sides an inherent rotation — that filaments that are dynamically 
cold should exhibit such a signal. 


One may draw an analogy with galaxies here: a disk galaxy 
viewed edge on shows a strong Azap (consider Azap as the 
Doppler shift from stars on the receding and approaching edge 
of the galaxy). The value (and significance as measured by the 
likelihood that it is random) of such a Azapg would decrease as 
the viewing angle of the galaxy is increased until the galaxy is 
seen face on. At this point the galaxy would appear “hot” - the 
motion of stars along the line of sight being random. Addition- 
ally, there is a degeneracy between morphology and shape in 
such situations: a very thin galaxy viewed with a large inclina- 
tion angle appears ellipsoidal. Since not all galaxies are rota- 
tionally supported, some galaxies that appear ellipsoidal will be 
inclined rotating thin disks while others will be dispersion sup- 
ported ellipsoids with little inherent rotation. It is this combina- 
tion of inclination angle and inherent rotation which explains the 
spectrum of Z-ms/Azap and significance seen in both the galaxy 
distribution as well as the spin of filaments. 


We note that it is the unique combination of small z,;m; and 
large Azap which returns a strong statistically significant signal. 
In other words there are filaments with large zrms and small Azap 
which clearly do not show any signal. This work does not pre- 
dict that every single filament in the universe is rotating, rather 
that there are subs-samples - intimately connected to the view- 
ing angle end point mass - which show a clear signal consistent 
with rotation. This is the main finding of this work. Again, we 
note that it is not trivial that filaments with low values of Zrms 
exhibit large values of Azag. This is only expected if filaments 
rotate. 


According to the zel’dovich approximation 5 the motion of 
galaxy’s is accompanied by larger scale matter flows #67) The 
mass and velocity flows are usually transported in a “hierarchi- 
cal” way: at first, matter collapses along the principle axis of 
compression forming great cosmic walls. Matter then flows in 
the wall plane along the direction of the intermediate axis of 
compression to form filaments. The final stage is the full 3 
dimensional collapse of an aspherical anisotrpoic density per- 
turbation wherein matter collapses and flows along the filament 
axes to form clusters. Under such a model of mass flow, one 
may expect that filament spin is formed at the second stage from 
the compression along the intermediate axis of compression. As 
such mass streams in the sheet towards the filament from both 
sides and could explain the observed signal. In other words the 
signal we measure is both consistent with a shear caused by 


galaxy motion in a sheet, or the spin of a filament embedded 
in that sheet. These two views are not contradictory with each 
other, since in the picture described above the spin of a filament 
is engendered by the shear of the sheet in which it is embedded. 

When the sample of filaments was divided into sub-samples 
based on the mass of the halos that sit at the ends of the fila- 
ments, those filaments with the most massive clusters have more 
significantly spinning signals than either the full sample or the 
filaments abutting the least massive halos. This gives a hint to- 
wards what may be driving the filament spin signal. It is also 
consistent with the idea that the spin is due to tidal fields influ- 
enced by the presence of massive aspherical density perturba- 
tions. The existence of a mass dependency is also strong evi- 
dence that what we have observed is not a systematic effect due 
to (e.g.) cosmological expansion. Any such systematic effects 
on the measurement of the redshift difference and its associa- 
tion with a velocity difference, would not manifest as a mass 
dependency. The mass dependency opens the door to a number 
of important questions as to the origin of this rotation or shear, 
questions that will be examined in future (theoretical) work such 
as: how is the measured effect engendered by the formation of 
clusters abutting filaments? Is there an angular momentum cor- 
relation between spinning filaments and nearby clusters? and 
perhaps most importantly: at which stage in structure forma- 
tion do filaments spin up, and how does this affect the galaxies 
within them? 

We know that from studies of the spatial distributions of 
smaller galaxies that inhabit the regions between large galax- 
ies, that they are accreted in a two fold process: first onto the 
filament connecting the two larger halos and then secondly an 
accretion along the line connecting the halos. In other words 
galaxies are first accreted from a direction perpendicular to the 
filament spine BSES] and then along the spine E8290] Assuming 
filaments can be approximated by a cylindrical potential such 
a flow would naturally give rise to an rotation about the fila- 
ment axis, followed possibly by helical accretion. This is con- 
sistent with the theoretical picture formed by #21 who’s analy- 
sis of numerical simulations indicate that the angular momen- 
tum of filaments is generated by those particles which collapse 
along the filament’s potential. Taken together, this paper and 221 
demonstrate that angular momentum can be generated on un- 
precedented scales, opening the door to a new understanding of 
cosmic spin. 
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Figure 1 — The statistical significance of Azag being consistent with random is shown as a function of the filament dynamical 
“temperature”, Zrms/Azag, in which Azap is the redshift difference of galaxies between the approaching and receding regions of 
each filament. The higher this quantity, the more unlikely it is that Azap is a random occurrence. Each filament is colored codes 
by the number of galaxies it contains, with yellow being rich filaments and purple being poor filaments as designated by the color 
bar. At a given value of Zms /Azap, the Azag seen in richer filaments is more statistically inconsistent with random. ‘Cold’ or ‘hot’ 
filaments are separated by Zrms/Azap = 1. 
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Figure 2 — The cumulative distribution of Azag , the redshift difference of galaxies in the approaching and receding regions of 
each filament. Left panel: all filaments irrespective of viewing angle. Middle panel: filaments whose axis subtends an angle 
cos@ < 0.2 with the line of sight. Right panel: filaments whose axis subtends an angle cos@ < 0.5 with the line of sight and 
which are dynamically cold, namely zms/Azag < 1. The red solid line shows the distribution of observed filaments, while the 
10,000 grey lines indicate the expected distribution from randomized redshifts. The median value of these 10,000 random samples 
is shown as the black solid line. The upper panels measure, as a function of Azap , the distance, in units of the standard deviation 
of the randomized distribution, between the measured curve and the mean of the random distributions. The upper x-axis displays 
the rotation velocity of the filament in km/s calculated as u = c x AZAB. 
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Figure 3 — The stacked rotation signal of filaments. Galaxies’ position rgaj long the filament are scaled by the length r; of the 
filament. dga is the distance of galaxies to the filament axis. Region A (defined as the region with greater mean redshift) is shown 
in the upper part of each plot, while region B is shown in the lower part. From top to bottom, each row shown the stacked rotation 
signal for all filaments, filaments whose spine subtend an angle cos(@) < 0.5 with the line of sight and filament with cos(@) < 0.2 
and which have Z;ms/Azap < 1. These three samples of filaments are subdivided according the group mass they are point to. 
Middle and right columns show the stacked rotation signal for filaments with the smallest and largest 10 percentile filament end 
point mass. The redshift difference is coded in color bar. 
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Figure 4 — The rotation curve of filaments as a function of the distance between galaxies to filament spine. The rotation speed 
is calculate by c x Az;, where Az; is the redshift difference of galaxies at given distance with respect to the mean redshift of all 
galaxies in the filament. The distance of galaxies in the region A (or receding region, defined as the region with greater mean 
redshift) is marked as positive values, while the distance of galaxies in the region B (or approaching region, defined as the region 
with smaller mean redshift) is marked as negative values. Error bars represent the standard deviation about the mean. 


METHODS 


The method presented here is fairly straightforward. We begin 
by using a publicly available catalog of filaments constructed from 
the SDSS DR12 galaxy survey Such a catalog is built using 
the Bisous algorithm 43 which identifies curvi-linear structures in the 
galaxy distributio 5355] Each filament is approximated by a cylin- 
der with a given axis that makes a viewing angle @ with the line of 
sight. Our hypothesis rests on the assumption that if filaments ro- 
tate then there should be a statistically significant component of a 
galaxy’s velocity that is perpendicular to the filament spine. Fur- 
thermore, this component should have opposite signs on either side 
of the filament spine: one receding and one approaching. In order 
to detect such a signal in principle we would only want to examine 
filaments whose spine is perpendicular to the line of sight such that 
a galaxy’s redshift (a proxy for the radial component of its veloc- 
ity) would coincide with its angular velocity. However restricting 
the filament sample this way dramatically reduces the sample size. 
Therefore we divide the filament catalog according to viewing an- 
gle into three bins: all filaments irrespective of viewing angle, those 
inclined by more than 60 degrees (cos@ < 0.5) with respect to the 
line of sight, those dynamically cold (see blow) and viewing angle 
more than ~80 degrees (cos @ < 0.2). This last sample is chosen for 
the reason that as demonstrated in Figure[A]] filaments inclined by 
more than ~80 degrees are perpendicular enough that the measured 
value of Zms/Azag is likely a good approximation to the inherent 
one. A consequence of this is there are more cold filaments with 
zms/Azap < 1 than hot filaments with zms/Azap > 1 inclined by 
cos@ < 0.2. In the appendix section the effect of viewing angle is 
examined in more detail. We find that the filaments whose spines 
are closer to perpendicular to the line of sight have more significant 
rotation signals (see figure A1). 

Each filament axis delineates two regions on either side of the fil- 
ament spine. The Bisous filament’s have an inherent scale of around 
1 Mpc. However all galaxies which are within a 2 Mpc distance from 
the filament axis are considered to be in one of these two regions. 
Galaxies further afield are not considered. We note that most of the 
galaxies that are within 2 Mpc are within 1Mpc and that including 
galaxies further away from the filament spine will “dilute” the signal 
as the regions further from the filament axis will, by construction, 
include random interlopers. 

The mean redshift of galaxies in both regions is calculated. The 
region with the greater mean redshift is arbitrarily called region A, 
the region with smaller mean redshift is called region B. The redshift 
difference Azag =(Za)-(Zp) is computed. The relative speed is then 
simply u = cAzap , where c is the speed of light. Next we exam- 
ine if the filament is dynamically hot or cold by using the redshift 
z of all galaxies in each filament to compute the rms namely Zms. 
This is then compared with Azap . Dynamically hot filaments have 
Zms >Azap while cold filaments have rms values less than the red- 
shift difference. Such a definition precludes the ability to search for a 
signal consistent with rotation in hot filaments since it is by definition 
flawed to ascribe a physical significance to a Azap if it is less than 
Zrms- Thus we do not expect to be able to find a discernible rotation 
signal in dynamically hot filaments. 

The median value of z;ms/Azap, (as well the standard deviation) 
is shown in Figure[AT]as a function of the angle subtended between 
the filament spine and the line of sight. Two important points can 
be observed by this plot: (1) the filament’s dynamical temperature 
(Zrms/AzZap) increases as filaments become parallel to the line of 


'The Bisous algorithm is not astronomy specific and was developed to iden- 
tify roads in satellite imaging. 


sight. Because the actual value of zms/Azap for a given filament 
is independent of @, this fact demonstrates that the viewing angle is 
critical in measuring this quantity and that z>ms/Azap can indeed be 
used to measure rotation for filaments that are close to perpendicular 
to the line of sight. Furthermore such a relationship is natural if a 
significant component of a galaxy’s velocity is perpendicular to the 
filament spine and consistent with rotation or shear. (2) The standard 
deviation of z:ms/Azap increases as filaments become parallel to the 
line of sight. This expected if in addition to rotational motion about 
the filament spine, galaxies are also traveling along the filament axis 
- perhaps with helical motion. 

A lower limit on the number of galaxies in each region of A and 
B of each filament must be assumed. Without such a limit we would 
introduce pathological cases into the analysis, namely filaments with 
no or just one or two galaxies in region A or B. Therefore we only 
consider filaments with at least 3 galaxies in each region - it is diffi- 
cult to justify the computation of an rms for anything less than this. 
The three galaxy per region A and B is a conservative choice and is 
the least arbitrary limit that can be applied. Despite this choice, the 
signal is examined as a function of galaxy number in the appendix 
section. In sum: increasing the minimum number of galaxies in each 
region decreases the sample size and increases the significance of 
the measured signal.The final fiducial sample consists of 17,181 fil- 
aments, and 213,625 galaxies in total. 

A Statistical significance to each measured value of Azap must 
be obtained to ensure that the measured signal is not simply the result 
of stochasticity. In other words we wish to ask the question “What 
is the chance that the observed distribution or Azap could be repro- 
duced given a random distribution of redshifts (in a given filament)?” 
In order to answer this question, the following test is performed. The 
redshifts of all galaxies in a given filament are randomly shuffled, 
keeping the position fixed. This is performed 10,000 times per fil- 
ament}? For each of these trials, the procedure described above is 
carried out: namely two regions are defined, their mean redshifts are 
measured and Az'andom is computed. The “real” Azap can then be 
compared to the distribution of random Aen on a filament by fil- 
ament basis. The chance that a measured Azag is consistent with 
random can then be quantified by where (how many standard devi- 


ations from the mean) the measured Azag is, in the distribution of 
Aztandom 
AB - 


A The effect of filament viewing angle and galaxy 
number 

In this section we examine how the signal changes based on two 
arbitrary choices we are forced to take in this work: the orientation 
of filaments with respect to the line of sight and the lower limit for 
the number of galaxies in each region to be considered. 

In order to examine how the orientation of filaments changes the 
significance of the signal we wish to measure, we perform the fol- 
lowing test. Let @ be the angle subtended between the filament spine 
and the line of sight. We compare filaments where cos @ < 0.05 since 
these are nearly perpendicular to the line of sight with those where 
cos @ > 0.95 corresponding to filaments nearly parallel to the line of 
sight. These are shown by the blue and orange histograms (respec- 
tively) in Figure Predictably we find that the filaments whose 
spines are closer to perpendicular to the line of sight have more sig- 
nificant rotation signals. Since the likelihood that a filament rotates 
does not depend on the orientation of the filament, we may ascribe 


*Tn the poorest filaments with only 6 galaxies (say), there are less than 10,000 
unique permutations. However performing this test 10,000 times, although re- 
dundant, will not affect the determination of the significance. 
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the decrease in significance to projection effects. 

Figure [A2] left, middle, and right panels consider the effect of 
changing the number of galaxies in each filament region A and B for 
3, 4, and 5 galaxies respectively. Note that because each filament 
has two regions so this corresponds to a minimum of 6, 8 and 10 
galaxies in the filament. We note that changing Ngai has little effect 
on the distribution of the estimated significance. Hence the choice of 
Ngai > 3 used throughout the paper is not biasing the result towards 
or away from significant detections. 


Data Availability The galaxy data used in this paper is drawn from the 


publicly available SDSS DR12 (https://www.sdss.org/dri2/) and 
can be found here: The Bisous filament cata- 


logue can be found here: |http: //cosmodb.to.ee 


Code availability The codes used in this study are available from the cor- 
responding authors upon reasonable request. 
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Figure A1 — The median dynamical “temperature”, Z:ms/Azap, of filaments as a function of the inclination angle @ between the 
filament spine and the line of sight. zms is the rms of galaxy redshift, and Azag is the mean redshift difference of galaxies on 
either side of the filament spine. That this quantity increases as the filaments become parallel to the line of sight is consistent with 
galaxies moving predominantly around the filament spine. Error bars represent the standard deviation about the median. 
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Figure A2 — The distribution of significance (with respect to the random sample) for filaments that are inclined nearly parallel 
(blue) and nearly perpendicular (orange) to the line of sight. In the left, middle right panel we show this for thre choices of 


Neal = 3, 4, 5. 
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